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PACS numbers: 74.25.Nf, 74.20. Rp 

We measure changes in the penetration depth A of the 
T c « 6 K superconductor LaFePO. In the process scanning 
SQUID susceptometry is demonstrated as a technique for 
accurately measuring local temperature-dependent changes 
in A, making it ideal for studying early or difficult-to-grow 
materials. A of LaFePO is found to vary linearly with tem- 
perature from 0.36 to ~2 K, with a slope of 143±15 A/K, 
suggesting line nodes in the superconducting order param- 
eter. The linear dependence up to ~ T c /3 is similar to the 
cuprate superconductors, indicating well-developed nodes. 

Research on the iron pnictide superconductors has 
been intense over the past year. Most attention has fo- 
cused on arsenic-based materials, which have the high- 
est transition temperatures, but which only supercon- 
duct at ambient pressure when doped, resulting in in- 
trinsic disorder. Superconductivity in LaFePO was an- 
nounced in 2006 pQ. Most likely it is the fully stoichio- 
metric compound that superconducts, with T c w 6 K; 
very clean residual resistivities below 0.1 m^-cm have 
been obtained [2]. How similar LaFePO will prove to 
be to the higher- T c arsenide compounds is not clear; al- 
though LaFePO does not show the magnetic order found 
in the arsenide compounds, its electronic structure has 
been found to be very similar [3]. 

The temperature dependence of the penetration depth 
provides information on the superconducting order pa- 
rameter (OP). OPs with line nodes are known to result 
in a T-linear dependence of AA(T) = A(T) — A(0) at 
low T [4 . Scattering modifies this dependence to T 2 [5] . 
A fully-gapped OP results in an exponential dependence 
AA oc T- 1 / 2 exp(-T /T) [6 . For the iron pnictide su- 
perconductors, proposed OPs include nodal and nodeless 

s [aiannniin], s+d nanu andd ^mcGani]. 

Most of these predictions are based on calculations in an 
unfolded, 1 Fe per unit cell Brillouin zone, which neglects 
the avoided crossings between the electron pockets in the 
true zone [3 j. These avoided crossings could significantly 
alter the nodal structure of the OP. It is also a possi- 
bility that different pnictide superconductors, although 
electronically similar, have different OPs [T7] . 

Radio- frequency tunnel diode resonator and microwave 



cavity perturbation measurements on iron arsenide su- 
perconductors have shown both power-law and exponen- 
tial temperature dependences of A A. Power-law depen- 
dence has been found in Ba(Fei_ x Co :E )2As2 [H], with 
the exponent n varying between 2.0 and 2.6 with dop- 
ing, n « 2 has been found in Bai_ x K :c Fe2As2 [19] , 
NdFeAsOo.gFo.i [20^ and LaFeAsOo. 9 F .i [20 . Expo- 
nential behavior has been found in Bai_ :E K :c Fe2As2 [21], 
PrFeAsOi-^ [22 and SmFeAsOo. 8 F .2 [23]. 

In LaFePO, nearly linear dependence of A A on T to 
below 150 mK has been reported by Fletcher et al [24] . 
using an RF tunnel diode circuit. However early LaFePO 
samples have had irregular shapes, which complicate RF 
and microwave measurements: to isolate A a 5 the mag- 
netic field of the excitation must be specifically oriented 
relative to the crystal axes, and at lower frequencies 
knowledge of the sample size is necessary to extract A A. 
Fletcher et al report these slopes dX/dT on three sam- 
ples: 412, 436 and 265 A/K (over 0.7 < T < 1.0 K). The 
magnitude of dX/dT constrains the number and opening 
angle of nodes, so confirmation with additional measure- 
ment is desirable. 

SQUID susceptometry has been demonstrated as a 
technique for observing superconducting transitions [25] , 
and has been used to determine the Pearl length A of thin 
superconducting films, for A ~10-100 /im [26 . We ex- 
tend this technique to measurement of nm-scale changes 
in local A with varying sample temperature. Our suscep- 
tometer is a niobium-based design [28] : its front end is 
shown in Fig. [I] The pick-up loop is part of a SQUID, and 
an excitation current (in this work, at 1071 Hz) is applied 
to the field coil. What is measured is the field coil-pick- 
up loop mutual inductance M. The susceptometer chip 
is polished to a point, aligned at an angle relative to the 
sample (in this work, « 16°), and mounted onto a 3-axis 
scanner. The Meissner response of superconducting sam- 
ples partially shields the field coil, so M decreases as the 
susceptometer approaches the sample. 

The schematic in Fig. [I] shows a model of the suscep- 
tometer. The field coil is taken as a wire loop of radius 
R at a height h above the superconductor surface, and 
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FIG. 1: Left: Field coil-pick-up loop mutual inductance M 
against height above the sample; ft CO ntact is the height at which 
the corner of the susceptometer chip contacts the sample. The 
line is a fit of eq. 1. Right: photograph of the front end of the 
susceptometer, and a schematic of the susceptometer-sample 
geometry. The dashed loop is the image field coil. 



the superconductor response field as an image coil placed 
a height 2/i e ff beneath the field coil, where the effective 
height h e fi = ft + A. The flux through the pick-up coil 
(radius a) is taken as the field at its center times its area. 
All coils are taken parallel to the surface, neglecting the 
alignment angle. These approximations give a conversion 
between M and ft e ff: 
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To measure changes in A the susceptometer is placed in 
contact with a flat afr-plane area of the sample, and the 
sample temperature T is varied. The contact is sufficient 
to overcome system vibration but weak enough to avoid 
excessive thermal coupling (the susceptometer is main- 
tained at « 0.3 K). The contact keeps ft constant, so 
changes to ft e ff are changes in X^: we are using the fact 
that, for ft ^> A, the response field of the superconductor 
is a function of ft e fr alone [27] . In this sense, the physical 
origin of eq. 1 is irrelevant as long as it accurately models 
the dependence of M on ft, which Fig. [I] shows to be the 
case. R and a are fitting parameters; they approximately 
match the actual dimensions of the susceptometer, but 
with precise values that vary with alignment angle and 
sample surface orientation; R and a are obtained sepa- 
rately for each sample. Crucial to this measurement, if 
the susceptometer is over a flat ab surface then the rel- 
evant penetration depth is A a 6 alone, even with nonzero 
alignment angle [27]. Due to the alignment angle, the 
minimum ft is ft CO ntact ~ 3 jam. 

What is the accuracy of measurement of AA? The fit to 
eq. 1 returns R and a consistent with a particular conver- 
sion constant, in /xm/V, between ft and applied voltage 
to the scanner, which is measured separately, in this work 
with ±5% accuracy. All AA quoted in this work have this 
±5% systematic uncertainty. Also, deviations from the 
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FIG. 2: AA(T) = A(T) - A(0) of Pb. A drift has been 
subtracted as described in the text. Inset: Open symbols: 
measurement of Gasparovic and McLean [29]. Filled: present 
data; the vertical error bars are the systematic ±5% error on 
all AA data in this Letter. 



fit give errors on A A up to 1.5%. At large A the assump- 
tion that the response field is a function of ft + A breaks 
down; numerical simulation shows that, at ft = 3 /im and 
A(0) = 5000 A, this assumption leads A A to be under- 
estimated by 1% at AA = 5000 A and 4% at 10,000 A. 
Thermal gradients from the susceptometer-sample con- 
tact have minimal effect: control tests on sapphire show 
that the change in ft e ff attributable to these gradients is 
no more than - 20 A for T varying between 1 and 8 K. 
By tracking T c of LaFePO, we determine that contact 
locally cools the sample by only ~40 mK at T = 6 K. 

As a test we measure the penetration depth of a lump 
of industrial-grade lead; the results and comparison with 
an earlier measurement are shown in Fig. [2] A ~ 100 A- 
scale downward drift of ft e ff , due to the sensor gradually 
pressing a dent into the soft lead surface, is subtracted 
from our data. The drift rate is T-independent and was 
measured separately from the data in Fig. [2] so the flat- 
ness of AA at low T is real. 

Fig. [3] shows the main result of this work: AA a 5 
vs. T for two LaFePO crystals (at the points indi- 
cated in Fig. [4je) and (f)). For both data sets AA was 
recorded over multiple temperature sweeps, both warm- 
ing and cooling, and found to follow the same path. A 
is seen to vary nearly linearly with temperature. Fitting 
AA = A + BT n over 0.7 < T < 1.6 K, from top to bot- 
tom n=1.22(4), 1.13(10) and 0.97(5) are obtained for the 
three curves in Fig. |3|a). 

Photographs of the two LaFePO specimens are shown 
in Fig. [4j An example of a susceptibility scan (a scan 
of the spatial variation in M) is shown in Fig. Qc). Be- 
cause M varies strongly with ft, features in individual 
scans mainly reflect surface topography. More useful is 
comparison of scans at different T: e.g. Fig. |4jcl) shows 
a map of ft e ff(3 K) — ft e fj(0.4 K) on sample #2, which 




FIG. 3: Top: A A of two LaFePO specimens, at the points in- 
dicated in Fig.fi] The lines were fit over 0.7 < T < 1.6 K. Bot- 
tom: black lines are possible superfluid densities for LaFePO 
sample #1, point 2, with different A(0). Shaded area: super- 
fluid density of YBa2Cus06.99 (l/(A a Ab)), from [30] and [3T] ; 
the width of the shaded area reflects uncertainty in A(0). 



reveals two useful facts: (1) Where the sample surface is 
not flat A c mixes in strongly and Ah e f? is large; one needs 
to be at least ~ 10 jam from an edge to measure A a 5. (2) 
Where the sample is flat, and Ah e ff = AA a 5, AA a 5 is 
homogeneous to within ~ 5%; areas of moderately in- 
creased Ah e ^ are areas where the surface is pitted. 

Maps of local T c , shown in Fig.[4^e) and (f), are made 
by performing susceptibility scans at various T and ex- 
tracting A/i e ff(T). Most areas of the samples show weak 
tails of superfluid density persisting a few 0.1 K beyond 
the dominant local T c (and in places beyond 7 K), which 
give uncertainty to estimates of the dominant T c . Our 
criterion for determining local T c is based on superfluid 
density: the A(0) = 4500 A superfluid density curve in 
Fig. |3|b) is taken as a reference, and in the scans the 
local A/i e ff is taken as A A. The local A(0) (which varies 
with topography) and T c are varied to obtain the best fit 
to the reference. Varying the reference A a 5(0) by 1000 A 
varies the calculated T c 's by ~ 0.1 K. 

A(0) could in principle be extracted from the geome- 
try of the susceptometer and its contact with the sample 
surfaces; however the uncertainties are large. From the 
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FIG. 4: (a,b) LaFePO specimens #1 and #2. (c) Suscep- 
tibility scan of #1 at T — 0.4 K. Over the superconductor 
M is reduced from its vacuum level by the Meissner response, 
(d) Change in /i e ff = h + \ between 0.4 and 3 K over specimen 
#2. (e,f) Maps of local T c , over the same areas as in (c) and 
(d). The crosses indicate the points where AA(T) data were 
collected. 

variation of M with surface orientation and SEM im- 
ages of the susceptometer a plausible contact point on 
the susceptometer can be identified, and comparison with 
the lead specimen indicates that A a 5(0) of LaFePO likely 
falls in the range 3500-5500 A. 

At the five points on sample #1 indicated in Fig.[4^e), 
dX/dT over the linear portion of A A is 146, 139, 136, 150 
and 205 A/K, and at the single measurement point on 
sample #2, 142 A/K. The 205 A/K measurement was 
at a point with significant topography and can be ex- 
cluded. Taking into account the 5% and 1.5% uncertain- 
ties, dX/dT is 143±15 A/K. 

The superfluid densities ps = 1/A 2 of LaFePO and 
YBa2Cu3 O6.99 are compared in Fig. [3] The linear por- 
tion of ps persists to a similar fraction of T c in both mate- 
rials, indicating that the nodes in LaFePO must be well- 
formed, as in YBa2Cu3 06.99 — the magnitude of the gap 
on either side of the nodes must be similar. In contrast, 
accidental nodes in nodal s orders may result in very 
asymmetric + and — lobes [32]. Also apparent in Fig. [3] 
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FIG. 5: (a) Photograph and (b) map of local T c of a T c « 
2 K specimen. The relative error on local T c 's is ~ 0.1 K. 
(c) AA(T) at the three points indicated in (b). Because the 
sample surface is not flat, these measured A A include some 
fraction of A c . 



is that ps of LaFePO rises very sharply on cooling just 
below T c . If the pairing is mediated by magnetic fluctu- 
ations then such a sharp rise may result from a gapping- 
out of low- frequency, pair-breaking fluctuations [33] . 

An intriguing possibility of highly asymmetric nodal s 
orders is that scattering might lift the nodes altogether, 
resulting in exponential rather than T 2 dependence of 
AA at low T [32 . In this work we also studied a T c w 
2 K LaFePO specimen. The reason for the anomalously 
low T c is unclear; electron probe microanalysis shows no 
impurities (to the \% level), and no anomalies in the La, 
Fe and P concentrations (to within \->\-> and 2%). 

Fig. [5] a shows the 2 K specimen. Compared with the 
6 K samples T c varies more widely both on large and 
small length scales: at each point studied strong tails 
of superfluid density extend well above the dominant lo- 
cal T c . A A versus temperature was recorded at three 
locations. Fitting to a power law over T < 0.8 K, expo- 
nents n=2.3±0.1, 2.0±0.1 and 1.6±0.1 are obtained. The 
dominant local T c 's at these three locations are 2.1±0.2, 
2.0dz0.1 and 2.5±0.1 K, respectively; deviations from 
n=2 may in part reflect variation in the local T c (or local 
T c distribution). On data up to 0.8 K, power law fits per- 
form better than exponential fits. Within our precision 
and temperature limits, AA(T) is consistent with dirty 
nodal superconductivity, and with many of the measure- 
ments on As-based materials. 

In conclusion, we have observed a linear temperature 
dependence of AA a b(T) below ~ T c /3 in LaFePO and 



accurately measured its slope, 143±15 A/K, using a local 
technique. The large temperature range of linear \ a b(T) 
indicates well- formed nodes. 
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